Differential gene expression and miRNA regulatory network in coronary slow flow.

Coronary slow flow (CSF) is characterized by slow progression of coronary angiography without epicardial stenosis. The aim of this study was to explore the potential biomarkers and regulatory mechanism for CSF. Peripheral blood mononuclear cells from 3 cases of CSF and 3 healthy controls were collected for high-throughput sequencing of mRNA and miRNA, respectively. The differentially expressed mRNAs (DE-mRNAs) and miRNAs (DE-miRNAs) was identified. A total of 117 DE-mRNAs and 32 DE-miRNAs were obtained and they were mainly enriched in immune and inflammatory responses. Twenty-six DE-mRNAs were the predicted target genes for miRNAs by RAID, and then the regulatory network of 15 miRNAs were constructed. In addition, through the PPI network, we identified the three genes (FPR1, FPR2 and CXCR4) with larger degrees as hub genes. Among them, FPR1 was regulated by hsa-miR-342-3p, hsa-let-7c-5p and hsa-miR-197-3p and participated in the immune response. Finally, we validated the differential expression of hub genes and key miRNAs between 20 CSF and 20 control. Moreover, we found that miR-342-3p has a targeted regulatory relationship with FPR1, and their expression is negatively correlated. Then we established a hypoxia/reoxygenation (H/R) HUVEC model and detected FPR1, cell proliferation and apoptosis. Transfection with miR-342-3p mimics can significantly promote the proliferation of HUVEC under H/R conditions. FPR1 were associated with CSF as a biomarker and may be regulated by miR-342-3p potential biomarkers.

Long noncoding RNAs and mRNAs profiling in ovary during laying and broodiness in Taihe Black-Bone Silky Fowls (Gallus gallus Domesticus Brisson).

BACKGROUND: Broodiness significantly impacts poultry egg production, particularly notable in specific breeds such as the black-boned Silky, characterized by pronounced broodiness. An understanding of the alterations in ovarian signaling is essential for elucidating the mechanisms that influence broodiness. However, comparative research on the characteristics of long non-coding RNAs (lncRNAs) in the ovaries of broody chickens (BC) and high egg-laying chickens (GC) remains scant. In this investigation, we employed RNA sequencing to assess the ovarian transcriptomes, which include both lncRNAs and mRNAs, in eight Taihe Black-Bone Silky Fowls (TBsf), categorized into broody and high egg-laying groups. This study aims to provide a clearer understanding of the genetic underpinnings associated with broodiness and egg production. RESULTS: We have identified a total of 16,444 mRNAs and 18,756 lncRNAs, of which 349 mRNAs and 651 lncRNAs exhibited significantly different expression (DE) between the BC and GC groups. Furthermore, we have identified the cis-regulated and trans-regulated target genes of differentially abundant lncRNA transcripts and have constructed an lncRNA-mRNA trans-regulated interaction network linked to ovarian follicle development. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) annotation analyses have revealed that DE mRNAs and the target genes of DE lncRNAs are associated with pathways including neuroactive ligand-receptor interaction, CCR6 chemokine receptor binding, G-protein coupled receptor binding, cytokine-cytokine receptor interaction, and ECM-receptor interaction. CONCLUSION: Our research presents a comprehensive compilation of lncRNAs and mRNAs linked to ovarian development. Additionally, it establishes a predictive interaction network involving differentially abundant lncRNAs and differentially expressed genes (DEGs) within TBsf. This significantly contributes to our understanding of the intricate interactions between lncRNAs and genes governing brooding behavior.

Gene Co-expression Network Analysis to Mine Genes Involved in Plant Secondary Metabolite Biosynthesis and Regulation.

Plants produce diverse specialized metabolites (SMs) that do not participate in plant growth and development but help them adapt to various environmental conditions. In addition to aiding in plant adaptation, different SMs serve as active ingredients for pharmaceutical and cosmetics products. However, despite their significant role in plant adaptation and industrial importance, the genes involved in the biosynthesis and regulation of many SMs remain largely unknown. This hinders deciphering the specific role of SMs in plant adaptation and limits their industrial utilization. Since many SMs pathway genes are expected to act in tight association with each other within a coexpression network, the network biology approach, such as weighted gene coexpression network analysis, could be used to identify the unknown genes. This chapter describes a workflow for constructing a gene coexpression network to identify genes that could be associated with the biosynthesis and regulation of SMs.

Trans-eQTL mapping in gene sets identifies network effects of genetic variants.

Nearly all trait-associated variants identified in genome-wide association studies (GWASs) are noncoding. The cis regulatory effects of these variants have been extensively characterized, but how they affect gene regulation in trans has been the subject of fewer studies because of the difficulty in detecting trans-expression quantitative loci (eQTLs). We developed trans-PCO for detecting trans effects of genetic variants on gene networks. Our simulations demonstrate that trans-PCO substantially outperforms existing trans-eQTL mapping methods. We applied trans-PCO to two gene expression datasets from whole blood, DGN (N = 913) and eQTLGen (N = 31,684), and identified 14,985 high-quality trans-eSNP-module pairs associated with 197 co-expression gene modules and biological processes. We performed colocalization analyses between GWAS loci of 46 complex traits and the trans-eQTLs. We demonstrated that the identified trans effects can help us understand how trait-associated variants affect gene regulatory networks and biological pathways.

Identification of genes related to ribosomal proteins in colorectal cancer: exploring their potential as biomarkers, prognostic indicators, and therapeutic targets.

BACKGROUND: Colorectal cancer (CRC) ranks as the third most commonly diagnosed cancer in both females and males, underscoring the need for the identification of effective biomarkers. METHODS AND RESULTS: We assessed the expression levels of ribosomal proteins (RPs) at both mRNA and protein levels. Subsequently, leveraging the STRING database, we constructed a protein-protein interaction network and identified hub genes. The co-expression network of differentially expressed genes associated with CRC and their target hub RPs was constructed using the weighted gene co-expression network analysis algorithm. Gene ontology and molecular signatures database were conducted to gain insights into the biological roles of genes associated with the identified module. To confirm the results, the expression level of the candidate genes in the CRC samples compared to the adjacent healthy was evaluated by the RT-qPCR method. Our findings indicated that the genes related to RPs were predominantly enriched in biological processes associated with Myc Targets, Oxidative Phosphorylation, and cell proliferation. Also, results demonstrated that elevated levels of GRWD1, MCM5, IMP4, and RABEPK that related to RPs were associated with poor prognostic outcomes for CRC patients. Notably, IMP4 and RABEPK exhibited higher diagnostic value. Moreover, the expression of IMP4 and RABEPK showed a significant association with drug resistance using cancer cell line encyclopedia and genomics of drug sensitivity in cancer databases. Also, the results showed that the expression level of IMP4 and RABEPK in cancerous samples was significantly higher compared to the adjacent healthy ones. CONCLUSION: The general results of this study have shown that many genes related to RPs are increased in cancer and could be associated with the death rate of patients. We also highlighted the therapeutic and prognostic potentials of RPs genes in CRC.

Group 3 medulloblastoma transcriptional networks collapse under domain specific EP300/CBP inhibition.

Chemical discovery efforts commonly target individual protein domains. Many proteins, including the EP300/CBP histone acetyltransferases (HATs), contain several targetable domains. EP300/CBP are critical gene-regulatory targets in cancer, with existing high potency inhibitors of either the catalytic HAT domain or protein-binding bromodomain (BRD). A domain-specific inhibitory approach to multidomain-containing proteins may identify exceptional-responding tumor types, thereby expanding a therapeutic index. Here, we discover that targeting EP300/CBP using the domain-specific inhibitors, A485 (HAT) or CCS1477 (BRD) have different effects in select tumor types. Group 3 medulloblastoma (G3MB) cells are especially sensitive to BRD, compared with HAT inhibition. Structurally, these effects are mediated by the difluorophenyl group in the catalytic core of CCS1477. Mechanistically, bromodomain inhibition causes rapid disruption of genetic dependency networks that are required for G3MB growth. These studies provide a domain-specific structural foundation for drug discovery efforts targeting EP300/CBP and identify a selective role for the EP300/CBP bromodomain in maintaining genetic dependency networks in G3MB.

Joint Bayesian estimation of cell dependence and gene associations in spatially resolved transcriptomic data.

Recent technologies such as spatial transcriptomics, enable the measurement of gene expressions at the single-cell level along with the spatial locations of these cells in the tissue. Spatial clustering of the cells provides valuable insights into the understanding of the functional organization of the tissue. However, most such clustering methods involve some dimension reduction that leads to a loss of the inherent dependency structure among genes at any spatial location in the tissue. This destroys valuable insights of gene co-expression patterns apart from possibly impacting spatial clustering performance. In spatial transcriptomics, the matrix-variate gene expression data, along with spatial coordinates of the single cells, provides information on both gene expression dependencies and cell spatial dependencies through its row and column covariances. In this work, we propose a joint Bayesian approach to simultaneously estimate these gene and spatial cell correlations. These estimates provide data summaries for downstream analyses. We illustrate our method with simulations and analysis of several real spatial transcriptomic datasets. Our work elucidates gene co-expression networks as well as clear spatial clustering patterns of the cells. Furthermore, our analysis reveals that downstream spatial-differential analysis may aid in the discovery of unknown cell types from known marker genes.

Single-cell insights: pioneering an integrated atlas of chromatin accessibility and transcriptomic landscapes in diabetic cardiomyopathy.

BACKGROUND: Diabetic cardiomyopathy (DCM) poses a growing health threat, elevating heart failure risk in diabetic individuals. Understanding DCM is crucial, with fibroblasts and endothelial cells playing pivotal roles in driving myocardial fibrosis and contributing to cardiac dysfunction. Advances in Multimodal single-cell profiling, such as scRNA-seq and scATAC-seq, provide deeper insights into DCM's unique cell states and molecular landscape for targeted therapeutic interventions. METHODS: Single-cell RNA and ATAC data from 10x Multiome libraries were processed using Cell Ranger ARC v2.0.1. Gene expression and ATAC data underwent Seurat and Signac filtration. Differential gene expression and accessible chromatin regions were identified. Transcription factor activity was estimated with chromVAR, and Cis-coaccessibility networks were calculated using Cicero. Coaccessibility connections were compared to the GeneHancer database. Gene Ontology analysis, biological process scoring, cell-cell communication analysis, and gene-motif correlation was performed to reveal intricate molecular changes. Immunofluorescent staining utilized various antibodies on paraffin-embedded tissues to verify the findings. RESULTS: This study integrated scRNA-seq and scATAC-seq data obtained from hearts of WT and DCM mice, elucidating molecular changes at the single-cell level throughout the diabetic cardiomyopathy progression. Robust and accurate clustering analysis of the integrated data revealed altered cell proportions, showcasing decreased endothelial cells and macrophages, coupled with increased fibroblasts and myocardial cells in the DCM group, indicating enhanced fibrosis and endothelial damage. Chromatin accessibility analysis unveiled unique patterns in cell types, with heightened transcriptional activity in myocardial cells. Subpopulation analysis highlighted distinct changes in cardiomyocytes and fibroblasts, emphasizing pathways related to fatty acid metabolism and cardiac contraction. Fibroblast-centered communication analysis identified interactions with endothelial cells, implicating VEGF receptors. Endothelial cell subpopulations exhibited altered gene expressions, emphasizing contraction and growth-related pathways. Candidate regulators, including Tcf21, Arnt, Stat5a, and Stat5b, were identified, suggesting their pivotal roles in DCM development. Immunofluorescence staining validated marker genes of cell subpopulations, confirming PDK4, PPARγ and Tpm1 as markers for metabolic pattern-altered cardiomyocytes, activated fibroblasts and endothelial cells with compromised proliferation. CONCLUSION: Our integrated scRNA-seq and scATAC-seq analysis unveils intricate cell states and molecular alterations in diabetic cardiomyopathy. Identified cell type-specific changes, transcription factors, and marker genes offer valuable insights. The study sheds light on potential therapeutic targets for DCM.

Identification of key mRNAs and signaling pathways in obsessive compulsive disorder based on weighted gene co-expression network analysis and cytoHubba plugin.

PURPOSE: Obsessive-compulsive disorder (OCD) is a complex psychiatric disorder. Genetic and broad environmental factors are common risk factors for OCD. The purpose of this study is to explore the molecular mechanism of OCD and to find new molecular targets for the diagnosis and management of OCD. METHODS: All data were downloaded from public dataset. Key modules and candidate key mRNAs were identified based on weighted gene co-expression network analysis (WGCNA). The "limma" R package was used for differential expression analysis of mRNAs. Subsequently, functional enrichment analysis of differentially expressed mRNAs (DEmRNAs) was also carried out. In addition, a diagnostic model was constructed. Finally, the infiltration level of immune cells in OCD and its correlation with multicentric key DEmRNAs were analyzed. RESULTS: Green and red modules were selected as the hub modules. A total of 447 mRNAs were considered candidate key mRNAs according to GS > 0.2 and MM > 0.3. A total of 26 DEmRNAs in the same direction were identified in the GSE60190 and GSE78104 datasets. A total of 26 DEmRNAs were intersected with candidate key mRNAs in WGCNA to obtain 10 intersection DEmRNAs (HSPB1, ITPK1, CBX7, PPP1R10, TAOK1, PISD, MKNK2, RWDD1, PPA1, and RELN). However, only four DEmRNAs (HSPB1, TAOK1, MKNK2, and PPA1) predicted related drugs. Subsequently, receiver operating characteristic analysis shows that the diagnostic model has high diagnostic value. Moreover, six multicentric key DEmRNAs (SNRPF, SNRNP70, PRPF8, NOP56, EPRS, and CCT2) were screened by UpSet package. Finally, six multicentric key DEmRNAs were found to be associated with immune cells. CONCLUSION: The key molecules obtained in this study lay a foundation for further research on the molecular mechanism of OCD.

Identification and bioinformatics analysis of lncRNAs in serum of patients with ankylosing spondylitis.

OBJECTIVES: The aim of this study was to explore the long non-coding RNA (lncRNA) expression profiles in serum of patients with ankylosing spondylitis (AS). The role of these lncRNAs in this complex autoimmune situation needs to be evaluated. METHODS: We used high-throughput whole-transcriptome sequencing to generate sequencing data from three patients with AS and three normal controls (NC). Then, we performed bioinformatics analyses to identify the functional and biological processes associated with differentially expressed lncRNAs (DElncRNAs). We confirmed the validity of our RNA-seq data by assessing the expression of eight lncRNAs via quantitative reverse transcription polymerase chain reaction (qRT-PCR) in 20 AS and 20 NC samples. We measured the correlation between the expression levels of lncRNAs and patient clinical index values using the Spearman correlation test. RESULTS: We identified 72 significantly upregulated and 73 significantly downregulated lncRNAs in AS patients compared to NC. qRT-PCR was performed to validate the expression of selected DElncRNAs; the results demonstrated that the expression levels of MALAT1:24, NBR2:9, lnc-DLK1-35:13, lnc-LARP1-1:1, lnc-AIPL1-1:7, and lnc-SLC12A7-1:16 were consistent with the sequencing analysis results. Enrichment analysis showed that DElncRNAs mainly participated in the immune and inflammatory responses pathways, such as regulation of protein ubiquitination, major histocompatibility complex class I-mediated antigen processing and presentation, MAPkinase activation, and interleukin-17 signaling pathways. In addition, a competing endogenous RNA network was constructed to determine the interaction among the lncRNAs, microRNAs, and mRNAs based on the confirmed lncRNAs (MALAT1:24 and NBR2:9). We further found the expression of MALAT1:24 and NBR2:9 to be positively correlated with disease severity. CONCLUSION: Taken together, our study presents a comprehensive overview of lncRNAs in the serum of AS patients, thereby contributing novel perspectives on the underlying pathogenic mechanisms of this condition. In addition, our study predicted MALAT1 has the potential to be deeply involved in the pathogenesis of AS.

Robust identification of interactions between heat-stress responsive genes in the chicken brain using Bayesian networks and augmented expression data.

Bayesian networks represent a useful tool to explore interactions within biological systems. The aims of this study were to identify a reduced number of genes associated with a stress condition in chickens (Gallus gallus) and to unravel their interactions by implementing a Bayesian network approach. Initially, one publicly available dataset (3 control vs. 3 heat-stressed chickens) was used to identify the stress signal, represented by 25 differentially expressed genes (DEGs). The dataset was augmented by looking for the 25 DEGs in other four publicly available databases. Bayesian network algorithms were used to discover the informative relationships between the DEGs. Only ten out of the 25 DEGs displayed interactions. Four of them were Heat Shock Proteins that could be playing a key role, especially under stress conditions, where maintaining the correct functioning of the cell machinery might be crucial. One of the DEGs is an open reading frame whose function is yet unknown, highlighting the power of Bayesian networks in knowledge discovery. Identifying an initial stress signal, augmenting it by combining other databases, and finally learning the structure of Bayesian networks allowed us to find genes closely related to stress, with the possibility of further exploring the system in future studies.

Analyses of lncRNA and mRNA profiles in recurrent atrial fibrillation after catheter ablation.

BACKGROUND: Atrial fibrillation (AF) is the most common cardiac arrhythmia worldwide. Catheter ablation has become a crucial treatment for AF. However, there is a possibility of atrial fibrillation recurrence after catheter ablation. Our study sought to elucidate the role of lncRNA‒mRNA regulatory networks in late AF recurrence after catheter ablation. METHODS: We conducted RNA sequencing to profile the transcriptomes of 5 samples from the presence of recurrence after AF ablation (P-RAF) and 5 samples from the absence of recurrence after AF ablation (A-RAF). Differentially expressed genes (DEGs) and long noncoding RNAs (DE-lncRNAs) were analyzed using the DESeq2 R package. The functional correlations of the DEGs were assessed through Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses. A protein‒protein interaction (PPI) network was constructed using STRING and Cytoscape. We also established a lncRNA‒mRNA regulatory network between DE-lncRNAs and DEGs using BEDTools v2.1.2 software and the Pearson correlation coefficient method. To validate the high-throughput sequencing results of the hub genes, we conducted quantitative real-time polymerase chain reaction (qRT‒PCR) experiments. RESULTS: A total of 28,528 mRNAs and 42,333 lncRNAs were detected. A total of 96 DEGs and 203 DE-lncRNAs were identified between the two groups. GO analysis revealed that the DEGs were enriched in the biological processes (BPs) of "regulation of immune response" and "regulation of immune system process", the cellular components (CCs) of "extracellular matrix" and "cell‒cell junction", and the molecular functions (MFs) of "signaling adaptor activity" and "protein-macromolecule adaptor activity". According to the KEGG analysis, the DEGs were associated with the "PI3K-Akt signaling pathway" and "MAPK signaling pathway." Nine hub genes (MMP9, IGF2, FGFR1, HSPG2, GZMB, PEG10, GNLY, COL6A1, and KCNE3) were identified through the PPI network. lncRNA-TMEM51-AS1-201 was identified as a core regulator in the lncRNA‒mRNA regulatory network, suggesting its potential impact on the recurrence of AF after catheter ablation through the regulation of COL6A1, FGFR1, HSPG2, and IGF2. CONCLUSIONS: The recurrence of atrial fibrillation after catheter ablation may be associated with immune responses and fibrosis, with the extracellular matrix playing a crucial role. TMEM51-AS1-201 has been identified as a potential key target for AF recurrence after catheter ablation.

Bioinformatics analysis of the potentially functional circRNA-miRNA-mRNA network in breast cancer.

Breast cancer (BC) is the most common cancer among women with high morbidity and mortality. Therefore, new research is still needed for biomarker detection. GSE101124 and GSE182471 datasets were obtained from the Gene Expression Omnibus (GEO) database to evaluate differentially expressed circular RNAs (circRNAs). The Cancer Genome Atlas (TCGA) and Molecular Taxonomy of Breast Cancer International Consortium (METABRIC) databases were used to identify the significantly dysregulated microRNAs (miRNAs) and genes considering the Prediction Analysis of Microarray classification (PAM50). The circRNA-miRNA-mRNA relationship was investigated using the Cancer-Specific CircRNA, miRDB, miRTarBase, and miRWalk databases. The circRNA-miRNA-mRNA regulatory network was annotated using Gene Ontology (GO) analysis and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway database. The protein-protein interaction network was constructed by the STRING database and visualized by the Cytoscape tool. Then, raw miRNA data and genes were filtered using some selection criteria according to a specific expression level in PAM50 subgroups. A bottleneck method was utilized to obtain highly interacted hub genes using cytoHubba Cytoscape plugin. The Disease-Free Survival and Overall Survival analysis were performed for these hub genes, which are detected within the miRNA and circRNA axis in our study. We identified three circRNAs, three miRNAs, and eighteen candidate target genes that may play an important role in BC. In addition, it has been determined that these molecules can be useful in the classification of BC, especially in determining the basal-like breast cancer (BLBC) subtype. We conclude that hsa_circ_0000515/miR-486-5p/SDC1 axis may be an important biomarker candidate in distinguishing patients in the BLBC subgroup of BC.

Bayesian Regression Facilitates Quantitative Modeling of Cell Metabolism.

This paper presents Maud, a command-line application that implements Bayesian statistical inference for kinetic models of biochemical metabolic reaction networks. Maud takes into account quantitative information from omics experiments and background knowledge as well as structural information about kinetic mechanisms, regulatory interactions, and enzyme knockouts. Our paper reviews the existing options in this area, presents a case study illustrating how Maud can be used to analyze a metabolic network, and explains the biological, statistical, and computational design decisions underpinning Maud.

The dynamic landscape of enhancer-derived RNA during mouse early embryo development.

Enhancer-derived RNAs (eRNAs) play critical roles in diverse biological processes by facilitating their target gene expression. However, the abundance and function of eRNAs in early embryos are not clear. Here, we present a comprehensive eRNA atlas by systematically integrating publicly available datasets of mouse early embryos. We characterize the transcriptional and regulatory network of eRNAs and show that different embryo developmental stages have distinct eRNA expression and regulatory profiles. Paternal eRNAs are activated asymmetrically during zygotic genome activation (ZGA). Moreover, we identify an eRNA, MZGAe1, which plays an important function in regulating mouse ZGA and early embryo development. MZGAe1 knockdown leads to a developmental block from 2-cell embryo to blastocyst. We create an online data portal, M2ED2, to query and visualize eRNA expression and regulation. Our study thus provides a systematic landscape of eRNA and reveals the important role of eRNAs in regulating mouse early embryo development.

Development of a Novel Prognostic Model for Esophageal Squamous Cell Carcinoma: Insights into Immune Cell Interactions and Drug Sensitivity.

Esophageal squamous cell carcinoma (ESCC) presents a five-year survival rate below 20%, underscoring the need for improved prognostic markers. Our study analyzed ESCC-specific datasets to identify consistently differentially expressed genes. A Venn analysis followed by gene network interactions revealed 23 key genes, from which we built a prognostic model using the COX algorithm (p = 0.000245, 3-year AUC = 0.967). This model stratifies patients into risk groups, with high-risk individuals showing worse outcomes and lower chemotherapy sensitivity. Moreover, a link between risk scores and M2 macrophage infiltration, as well as significant correlations with immune checkpoint genes (e.g., SIGLEC15, PDCD1LG2, and HVCR2), was discovered. High-risk patients had lower Tumor Immune Dysfunction and Exclusion (TIDE) values, suggesting potential responsiveness to immune checkpoint blockade (ICB) therapy. Our efficient 23-gene prognostic model for ESCC indicates a dual utility in assessing prognosis and guiding therapeutic decisions, particularly in the context of ICB therapy for high-risk patients.

Tissue-specific atlas of trans-models for gene regulation elucidates complex regulation patterns.

BACKGROUND: Deciphering gene regulation is essential for understanding the underlying mechanisms of healthy and disease states. While the regulatory networks formed by transcription factors (TFs) and their target genes has been mostly studied with relation to cis effects such as in TF binding sites, we focused on trans effects of TFs on the expression of their transcribed genes and their potential mechanisms. RESULTS: We provide a comprehensive tissue-specific atlas, spanning 49 tissues of TF variations affecting gene expression through computational models considering two potential mechanisms, including combinatorial regulation by the expression of the TFs, and by genetic variants within the TF. We demonstrate that similarity between tissues based on our discovered genes corresponds to other types of tissue similarity. The genes affected by complex TF regulation, and their modelled TFs, were highly enriched for pharmacogenomic functions, while the TFs themselves were also enriched in several cancer and metabolic pathways. Additionally, genes that appear in multiple clusters are enriched for regulation of immune system while tissue clusters include cluster-specific genes that are enriched for biological functions and diseases previously associated with the tissues forming the cluster. Finally, our atlas exposes multilevel regulation across multiple tissues, where TFs regulate other TFs through the two tested mechanisms. CONCLUSIONS: Our tissue-specific atlas provides hierarchical tissue-specific trans genetic regulations that can be further studied for association with human phenotypes.

High expression of CASP1 induces atherosclerosis.

Atherosclerosis is a chronic, progressive vascular disease. The relationship between CASP1 gene expression and atherosclerosis remains unclear. The atherosclerosis dataset GSE132651 and GSE202625 profiles were downloaded from gene expression omnibus. Differentially expressed genes (DEGs) were screened. The construction and analysis of protein-protein interaction network, functional enrichment analysis, gene set enrichment analysis, and Comparative Toxicogenomics Database analysis were performed. Gene expression heatmap was drawn. TargetScan was used to screen miRNAs that regulate central DEG. 47 DEGs were identified. According to gene ontology analysis, they were mainly enriched in the regulation of stimulus response, response to organic matter, extracellular region, extracellular region, and the same protein binding. Kyoto Encyclopedia of Gene and Genome analysis results showed that the target cells were mainly enriched in the PI3K-Akt signaling pathway, Ras signaling pathway, and PPAR signaling pathway. In the enrichment project of Metascape, vascular development, regulation of body fluid levels, and positive regulation of cell motility can be seen in the gene ontology enrichment project. Eleven core genes (CASP1, NLRP3, MRC1, IRS1, PPARG, APOE, IL13, FGF2, CCR2, ICAM1, HIF1A) were obtained. IRS1, PPARG, APOE, FGF2, CCR2, and HIF1A genes are identified as core genes. Gene expression heatmap showed that CASP1 was highly expressed in atherosclerosis samples and low expressed in normal samples. NLRP3, MRC1, IRS1, PPARG, APOE, IL13, FGF2, CCR2, ICAM1, HIF1A were low expressed in atherosclerosis samples. CTD analysis showed that 5 genes (CASP1, NLRP3, CCR2, ICAM1, HIF1A) were found to be associated with pneumonia, inflammation, cardiac enlargement, and tumor invasiveness. CASP1 gene is highly expressed in atherosclerosis. The higher the CASP1 gene, the worse the prognosis.

Weighted gene co-expression network analysis reveals key biomarkers and immune infiltration characteristics for bronchial epithelial cells from asthmatic patients.

BACKGROUND: Asthma ranks among the most prevalent non-communicable diseases worldwide. Previous studies have elucidated the significant role of the immune system in its pathophysiology. Nevertheless, the immune-related mechanisms underlying asthma are complex and still inadequately understood. Thus, our objective was to investigate novel key biomarkers and immune infiltration characteristics associated with asthma by employing integrated bioinformatics tools. METHODS: In this study, we conducted a weighted gene co-expression network analysis (WGCNA) to identify key modules and genes potentially implicated in asthma. Functional annotation of these key modules and genes was carried out through gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis. Additionally, we constructed a protein-protein interaction (PPI) network using the STRING database to identify 10 hub genes. Furthermore, we evaluated the relative proportion of immune cells in bronchial epithelial cell samples from 20 healthy individuals and 88 asthmatic patients using CIBERSORT. Finally, we validated the hub genes and explored their correlation with immune infiltration. RESULTS: Furthermore, 20 gene expression modules and 10 hub genes were identified herein. Among them, complement component 3 (C3), prostaglandin I2 receptor (PTGIR), parathyroid hormone-like hormone (PTHLH), and C-X3-C motif chemokine ligand 1 (CX3CL1) were closely correlated with the infiltration of immune cells. They may be novel candidate biomarkers or therapeutic targets for asthma. Furthermore, B cells memory, and plasma cells might play an important role in immune cell infiltration after asthma. CONCLUSIONS: C3, PTGIR, CX3CL1, and PTHLH have important clinical diagnostic values and are correlated with infiltration of multiple immune cell types in asthma. These hub genes, B cells memory, and plasma cells may become important biological targets for therapeutic asthma drug screening and drug design.

CDK1 and CCNA2 play important roles in oral squamous cell carcinoma.

Oral squamous cell carcinoma (OSCC) is a malignant tumor that occurs in oral cavity and is dominated by squamous cells. The relationship between CDK1, CCNA2, and OSCC is still unclear. The OSCC datasets GSE74530 and GSE85195 configuration files were downloaded from the Gene Expression Omnibus (GEO) database and were derived from platforms GPL570 and GPL6480. Differentially expressed genes (DEGs) were screened. The weighted gene co-expression network analysis, functional enrichment analysis, gene set enrichment analysis, construction and analysis of protein-protein interaction (PPI) network, Comparative Toxicogenomics Database analysis were performed. Gene expression heatmap was drawn. TargetScan was used to screen miRNAs that regulate central DEGs. A total of 1756 DEGs were identified. According to Gene Ontology (GO) analysis, they were predominantly enriched in processes related to organic acid catabolic metabolism, centromeric, and chromosomal region condensation, and oxidoreductase activity. In Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis, the DEGs were mainly concentrated in metabolic pathways, P53 signaling pathway, and PPAR signaling pathway. Weighted gene co-expression network analysis was performed with a soft-thresholding power set at 9, leading to the identification of 6 core genes (BUB1B, CCNB1, KIF20A, CCNA2, CDCA8, CDK1). The gene expression heatmap revealed that core genes (CDK1, CCNA2) were highly expressed in OSCC samples. Comparative Toxicogenomics Database analysis demonstrated associations between the 6 genes (BUB1B, CCNB1, KIF20A, CCNA2, CDCA8, CDK1) and oral tumors, precancerous lesions, inflammation, immune system disorders, and tongue tumors. The associated miRNAs for CDK1 gene were hsa-miR-203a-3p.2, while for CCNA2 gene, they were hsa-miR-6766-3p, hsa-miR-4782-3p, and hsa-miR-219a-5p. CDK1 and CCNA2 are highly expressed in OSCC. The higher the expression of CDK1 and CCNA2, the worse the prognosis.